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PROCEDURES FOR 

HEURISTIC SCHEDULING UNDER LIMITED RESOURCES 
IN ACTIVITY NETWORKS 

Abstract 

This is the documentation of a program that permits the analyst 
to choose any one of three priority rules to schedule activities under 
conditions of limited resources. The program is written in FORTRAN. 

We give the logic of the programs, a flowchart representation, and a 
printout of the instructions (upon request). Also, a solved sample 
problem is included. 

This program is available at a nominal fee from the Graduate 
Program in Operations Research, N. C. State University, P. 0. Box 5511, 
Raleigh, NC 27650. 







PROCEDURES FOR HEURISTIC SCHEDULING UNDER 


LIMITED RESOURCES IN ACTIVITY NETWORKS 

S. E. Elmaghraby and Z. M. Naman 

This package consists of two programs: 

1. SETUP: To setup the variables in a form readable 

by HSP (heuristic scheculing procedure) 

2. HSP: To schedule the activities prepared by SETUP. 

I. SETUP 

This program is to set up the variables needed for the schedul¬ 
ing program HSP, namely the network, activities 1 duration, resource 
consumption and resources availabilities. The following are its options 
GN: GN=1 The network is randomly generated 
GN=2 The network is read in 

GRAND: GRAND=1 The network parameters are randomly generated. 
GRAND=2 The network parameters are read in 
CONV: C0NV=1 Convert the network to A-on-N. 

C0NV=2 The network is A-on-N. 

Any combination of the above options could be used. For 
example. 


GN=1 , GRAND=2 , C0NV=1 

means: To generate the network randomly; 

The network parameters are read in; and 

The network is A-on-A, so convert it to A-on-N. 
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Input Sequence 

1. Options: GN, GRAND, CONV 

2. Number of nodes, number of arcs 

3. IF_ the network is read THEN 

IF it is A-on-A THEN give a list of pairs (Node i. 
Node j) to represent arcs in the network 
ELSE [The network is A-on-N] give the upper triangular 
matrix without diagonal (the adjacency matrix) 
which is a (0,1) matrix given row by row. 

ELSE Generate the network randomly; 

4. Read the number of resources 

5. I_F the network parameters are randomly generated THEN give 

the following information: 


1 . 

L3: 

Maximum 

desired 

activity duration 

2. 

N3: 

Maximum 

desired 

resource requirement 

3. 

F > 1: 

Resource allocation factor 


[To get the resource allocation we find the average 
resource used, then we multiply it by F. 

Average resource = Total resource/critical path.] 
ELSE [The parameters are read in] give the following 
information: 

1. Duration of each activity 

2. Resource requirements to he given as FOR each 
activity 







FOR each resource 


R1: Resource required 

Dl: Duration of this request 

HND 

Resource availability for each resource 
END 

Example: 1, 1, 1 [option card] [generate the 

network, generate the parame¬ 
ters, convert to A-on-N] 

20, 60 [network specification card] 

r20-nodes and 60-arcs] 

3 [number of resources card] 

[3-resources] 

15, 10, 1,5 [parameters generated informa¬ 
tion] [L3=15, N3=10, F=1,5] 

Intermediate Output : 

1. Provided or generated information 

2. Network as A-on-N (arc list) 

3. Number of activities 

4. Upper bound on project duration 

5. Lower bound on project duration 

6. Activity data 

7. Resource availability levels. 

Program Limitations: 

1. Maximum number of resources = 4 

2. Maximum number of activities = 99 

3. Maximum sum of durations - 350 

To increase the above variables, the program DIMENSION 
must be changed. 


A 
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4. To add the network generator, do the following: 

a. Replace the comment card which indicates the 
position of the network generator by calling 
card. 

b. Take out all the "c" at Column 1 in the next 
Do-loop. 

c. Add the arrays NE and NS in the DIMENSION 
statement. 

5. For the random number generator, please look at the 
note given in ’’NETWORK GENERATION” program descrip¬ 
tion, OR Report No. 179, June 1981, p. 4. 

II. HSP 

The project scheduling problem is to assign start times to all 
the activities in A. Each of these start times must respect the pre¬ 
cedence and resource constraints described above. Another assumption 
made in this study is that of non-preemptive scheduling: Once started, 
an activity must continue without interruption for its duration. As 
well as respecting all constraints, schedulers attempt to assign start 
times so that some objective is achieved. 

The objective of the scheduling algorithm presented below is to 
minimize the project duration. Define the set, A(i), to be the set of 
activities immediately succeeding i. By L we denote the set of activi¬ 
ties having no successors. Activity i is in L if A(i) is the empty 

set. Project duration is then obtained as the max (t(i) + d(i) - 1) . 

ieL 

Necessarily, the maximization of total resource utilization is accom¬ 
plished by minimizing project duration. Moreover, minimizing project 









duration also assures the minimization of delay penalties incurred 
per period beyond the (resource free) critical path length. Finally, 
the chosen objective is applicable to a variety of situations. For 
example, a capital intensive project, such as a large construction 
project, might require substantial financing to secure the necessary 
resources. Interest payments on the borrowed capital cut into the 
constructor*s profit until the project is completed and the resources 
are no longer needed. Minimizing project duration thus minimizes 
such interest costs. 

Heuristics Employed 

In the heuristic scheduling procedure, an adapted bi npacking 
rule is used to assign the activities to the resource units. The 
first fit decreasing (FFD) binpacking rule is one of the most popular 
scheduling rules. It has also been applied with some success to 
the problem of memory storage allocation by some computer scientists. 
The FFD rule specifies the packing of different sized objects into 
containers or bins. The objects are first ordered by decreasing size. 
The largest object is then assigned to the first container in which 
it will fit. The second largest object is then assigned to the first 
container in which it will fit, etc. In the project setting, the 
objects are the jobs and the containers are the resources. There are 
as many "containers’ 1 as there are different resources, with each 
job "occupying" all containers simultaneously. This is a 
multi-dimensional binpacking problem which we wish to avoid by 


u 



adopting a simple rule for "sizing" the jobs. We address this 
"sizing" problem next. 

For project scheduling, it is clear that the limited avail¬ 
ability of the various resources may cause an activity to be delayed. 
Such delays can result in project completion delay. A measure of 
the delay in project completion caused by delaying an activity 
should serve well as an activity’s "size." 

The original idea for measuring the project delay caused by 
delaying activity i was to sum the durations of all activities 
succeeding activity i in the precedence ordering plus the duration 
of activity i itself. The larger the priority number ("size") of 
activity i, p(i), the longer the expected project delay should be, 
if activity i is delayed. To reduce the computational burden, we 
adopted a slightly different priority number: 


d(i) 


if A(i) = t. 

d(i) + 

Y pa) - 

if A(i) + 0. 


jeA(i) 


These numbers can be assigned in one pass through the precedence 
ordering. However, the recursion allows the counting of some dura¬ 
tions more than once in determining the priority number for certain 


activities. 
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The two other priority rankings available in this package 
are: the longest remaining path (LRP) heuristic and the shortest 
imminent activity (SIA) heuristic. LRP gives a higher priority to 
an activity with a longer critical subpath (to project completion) 
than a competing activity. The SIA heuristic mimics the shortest 
processing time rule from job shop scheduling. Characteristic of 
both schemes, and a number of other heuristics, is the neglect of 
the total impact on project completion that delaying an activity can 
produce. The impact of delaying an activity can affect all its 
successors. The proposed priority ranking is an attempt to reflect 
the full impact of delaying an activity on the project duration by 
adopting a simple surrogate measure. 

To sum up, the three heuristic "sizing" priority rules are: 

CR1: The sum of durations of all succeeding 
activities 

CR2: The longest path to project termination (LRP) 

CR3: The shortest imminent (available) activity (SIA) 

The scheduling process begins by forming the set, CS, of all 
activities where predecessors have been scheduled prior to the current 
scheduling period. On the first iteration, CS will consist of all 
initial activities (i such that B(i) = (j)). In general CS will contain 
activities that are unordered with respect to precedence. The adopted 
FFD rule assigns activities of CS to available resource units. The 
activities are the items, ranked by decreasing "size" (priority number). 
The "size" of each bin in a particular scheduling period will be the 
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number of available units of that resource type. The activities are 
checked for fit starting with that activity whose priority number is 
the largest in CS. If activity i "fits", then activity i is scheduled 
to begin on period Tl. After updating resource availabilities, the 
next activity is checked for fit, etc., until all activities in CS 
have been considered for scheduling at Tl, where T1 is the current 
scheduling period. 

The algorithm is iterative in nature. CS and Tl are updated 
after each iteration. The process repeats until all activities are 
assigned. The updated scheduling period is determined by the completion 
times of the activities scheduled during the previous iteration. It is 
the earliest completion time of any of the activities scheduled in the 
previous iteration. 

Tl = Tl . . + min (M(i, Tl • d(i)), 
new old ^ v * old' v 


where d(i) is the duration of activity i, and 


= < 


1 if activity i started on Tl 


00 if i not started on Tl , 

old 


The intention in skipping over the intervening periods between Tl ^ 
and Tl new is to avoid wasting time attempting to schedule activities 
when few resources are available in the scheduling period. Upon 
completion of an activity, resource units are "released" and other 
activities can perhaps be undertaken. A flowchart of HSP logic is 
shown in Fig. 1. 



An Example 


/(2,3) 



(6,3) 


(1,3) 


(3,3) 



(d(i), c (i) ) 


(7,3) 


(4,3) 


(5,3) 


The figure above is a graphical representation of the precedence 
relationship among the seven activities of an example project. An arc, 
(i,j) implies that activity j cannot be started until activity i has 
been completed. For example, activities two, three, and four cannot 
be started until activity one has been completed. Assume that activity 
durations, d(i) - i, for i = 1,2, ..., 7. Also, one resource type is 
required in quantities of three units per day by each activity 
(c^(j) = 3, for j = 1,2, ..., d(i)) . Six units of the resource are 
available for each of the scheduling periods (r(t) = 6, for 
t = 1,2, ..., T where T is an upper bound on project duration). 

First, the priority numbers are assigned in one backward pass 
through the precedence network. Suppose we adopt priority rule CR1. 
Activities six and seven are the terminal activities of this project 
(i.e., L = {6,7}). Thus their priority numbers are just their durations 
(p(6) = d(6) = 6, p(7) = d(7) = 7). Priorities for the immediate pre¬ 
decessors of six and seven, activities two, three, and five, can now 
be computed. Their priorities are the sum of activity duration and 
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the priority numbers of the successors: p(2) = d(2) + P (6) = 2 + 6 = 8, 
p(3) = d(3) + p(6) =3 + 6 = 9, p(5) = d(S) + p(7) = 5 + 7 = 12. 
Activity four has priority p(4), given by the sum of its duration and 
the priority number of its successor, activity five. Thus p(4) = 
d(4) + p(5) = 4 + 12 = 16. Finally, activity one has as successors 
activities two, three, and four. Thus p(l) = d(l) + p(2) ♦ p(3) + 
p(4) =1+8+9+ 16 =34. Now that activity priorities have been 
computed, the scheduling process can begin. 

On the first scheduling iteration, CS will consist of activity 
one since it is the lone initial activity in the project. Three units 
of the resource are required by activity one, and six units are avail¬ 
able on day one. Thus activity one is scheduled for day one. 

Day two is the updated scheduling period on the second iteration 
since activity one is completed on day one. CS will consist of activi¬ 
ties two, three, and four, the immediate successors of activity one. 

The scheduling priority is four (p(4) = 16), three (p(3) = 9), and 
two (p(2) = 8). Activity four is first checked for fit. Three units 
of the resource are required on days two through five, and six units 
are available. Since activity one is completed by day two, activity 
four is scheduled to start on day two. Now only three resource units 
remain available on days two through five. Next, activity three is 
checked for fit. Three resource units are required by activity three 
on days two through four, and three units are available. Activity one 
is completed by day two, so activity three is scheduled to start on 
day two. Thus no units of the resource are available on days two, 
three, and four. Activity two is checked for fit next. But, with no 



resources available and three units required, activity two cannot be 
scheduled on day two. So the set CS and scheduling period are updated 
and we proceed to the third iteration. 

The new set of activities to be scanned will include activity 
two that was not scheduled on the previous iteration and activity five 
since its predecessor has now been scheduled. The new scheduling period 
is day five since, of the two activities scheduled on the previous 
iteration, activity three finishes sooner (on day four) releasing some 
resources. Activity five has a larger priority number than activity 
two. Hence activity five is considered for scheduling before activity 
two. However, activity four is still in progress on day five so activ¬ 
ity five is ineligible for scheduling on day five. Next, activity two 
is considered for scheduling on day five. Three units of the resource 
are required by activity two on days five and six. Three units are 
available on day five, and six units are available on day six. No 
precedence constraint is violated, so activity two is scheduled to 
begin on day five. With that assignment, no resource units remain 
available on day five and three units remain available on day six. 

Both activities in CS, at this iteration, have been checked. We pro¬ 
ceed to update CS and the scheduling period for the next iteration. 

The updated CS consists of activities five and six since the scheduling 
of activity two releases activity six for consideration. The new 
scheduling period will be day six since activity four is completed on 
day five. 

On the fourth iteration, activity five has priority and is 
scheduled on day six. Available resource units are updated and 






activity six is considered. Sufficient resources are available, and 
no precedence constraint is violated in starting activity six on day 
seven. Hence, that assignment is made. 

Finally, activity seven is considered on the fifth iteration. 
The updated scheduling period is day eleven (given by the completion 
of activity five). No resource or precedence constraints are violated 
by starting activity seven on day eleven so that assignment is made. 
The resource "skyline” or profile chart below illustrates the example 
project’s schedule. The project is completed on day seventeen. This 
duration coincides favorably with the project’s critical path length 
in the absence of resource constraints (seventeen days), and is there¬ 
fore optimal. 



The project scheduling algorithm described above is an analog 


of the heuristic algorithm for the ALB problem called the "Ranked 
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Positional Weight" heuristic, which was proposed by Hegelson and Birnie 
in 1961. An extension of this heuristic was suggested by Mansoor in 
1964. 


Advantages of Heuristic Scheduling 

A significant contribution of heuristic scheduling might be that 
of increasing the computational efficiency of optimal project schedul¬ 
ing procedures. One attempt might be to link the proposed scheduling 
algorithm (or any other good heuristic algorithm) with a branch and 
bound algorithm. For such a tandem algorithm, the size of the space 
to be searched will be smaller than the search space of a straight¬ 
forward branch and bound procedure. 

We have noted that a lower bound on project duration is the 
critical path length. Similarly, we obtain a lower bound on the start 
time of an activity by computing the earliest start schedule (ESS). 

The ESS neglects any resource constraints and starts an activity as 
soon as all precursory activities are completed: 


ESS(i)- 


r 


< 


max 

jeB(i) 


(ESS(j) + d(j)). 


if B(i) 4 4 , 


1 if B(i) = 4. 


If a tight upper bound on optimal start time is available for all 
activities, then the chance of finding an optimal schedule by a 
branch and bound scheme is greatly improved. A good heuristic 
algorithm will provide such a bound at a small price. 










The heuristic’s schedule yields a feasible project duration, T. 


The late start schedule, LSS^, analogous to the early start schedule 
described above, is computed: 


LSS T (i) 


min (LSS (j) - d(i) + 1), if A(i) 4 <(l, 
jeA(i) 


T - d(i) + 1, if A(i) = (j>. 


The optimal start time(s) of activity i is (are) now bounded: 

ESS(i) <_ T t (i) <_ LSS^,(i) . Activity i cannot be started before its 

predecessors are completed (ESS(i)) and an optimal schedule cannot 

start activity i after its late start time with project duration T 

since T < T. 
opt — 

Thus the search tree of the branch and bound procedure is pruned. 
The tightness of these bounds and hence the reduction of the search 
space is increased as the heuristic schedule is improved. Therefore, 
an efficient, dependable heuristic scheduling algorithm 
might play a crucial role in increasing the computational reliability 
of an analytical procedure for the project scheduling problem. 

The following information is needed to be passed from 

SETUP: 

1. Number of resources 

2. Number of activities 

3. Upper bound of planning horizon 

4. Network in A-on-A 

5. Availability of each resource 

6. Requirement of each resource/activity/time period 







7. Duration of each activity 

8. Index vector K 


9. Early start 

schedule 

Options 



P7: 

P7=l 

Do not print iteration information. 


P7*2 

Print iteration information. 

SCH: 

SCH=1 

Scheduling using CR1 


SCH=2 

Scheduling using CR1 and CR2 


SCH=3 

Scheduling using CR1 and CR3 


SCH=4 

Scheduling using CR1, CR2, and CR3 


For example, P7=2 
SCH=3 


imply: Print iteration information. 

Use CR1 and CR3 as priority criteria. 

Input Sequence 

p 7, SCH [on the same card] 

Output 

1. For each scheduling iteration (optional) 

a. Current scanned subset of activities and their 
priority numbers 

b. Activities delayed because of resource limited 
availability 

c. Activities delayed because of precedence 

d. Activities scheduled 

2. Project table 

For each activity Print: 

Duration, Priority number, HSP-start, Early-start 

i 

:! 

1 
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3. Resource utilization 
For each resource Print: 

a. Time, availability, utilization, leftover 

b. Draw the skyline for the resource utilization 
vs . time. 

4. Adjusted percent resource utilization, for all resources 

5. Lower bound on duration 

6. Completion time from HSP 

7. Upper bound of duration 
Limitation 

Please see the limitation given in SETUP. 
































(OPTION) To read the 
options for this RUN. 

Read all the necessary 
information depending on 
the option that is given. 

P = Number of nodes if 
A-on-N 

= Number of arcs if 
A-on-A 

K = An array to tell where 
each resource require¬ 
ment is started for 
each activity. 


OPTION _ 

READ, GN, GRAND, CONY _ 

Print the option used in this run 
READ, NARC, NODE _ 

r\ gn = 


NETGEN 


K = 1 


C0NV=1 


NUP = 2*NARC NUP = NODE* 


FOR I = 1 to NARC READ, V(J) (N0DE-l)/2 


V(K) = NE(I) 

V (K+l) = NS (I) 
K = K + 2 


J = 1, NUP J = 1, NUP 


READ, RI 
GRAND = 2 
READ, L3, N3, F 
-*• CONVERT 


READ, D(I), I 

= 1, P 

I = 1 to P 


FOR RO = 1, RI 


READ, RI, D1 


FOR DO = 1, D 

1 

R(DO+K(I),R0 

I=Ri_ 


_ READ, A(1,R0 

FOR I - 1, P 

r- 1 " 1 .. - - — 

B2 = B2 ♦ PCI) 

FOR RO = 1, RI _ 

FOR I = 2, B2 
A(I,RO) = A(1,R0) 
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2. (CONVERT) To convert 
A-on-A to A-on-N and 
create an array X which is 
the one dimensional repre¬ 
sentation of the upper 
triangle w/o diagonal 
matrix. 

X = 1 if activity 
I precedes J 
= 0 otherwise 

3. (GEN) To generate randomly 
the network parameters. 

D(1) - duration of each 
activity 

R - resource requirement 
for each activity at 
each time period for 
each resource. 

4. (CPL) To calculate the 
critical path of the network. 
The early start schedule 
could be produced by 

E(J) = V(J) - D(J) + 1; 

1 J <_ P where V is the 
vector generated by CPL. 


CONVERT 


CONV = 2 


FOR I = 1, p - 1 


1X1 = (I-l)*(P-I/2) - I 


FOR J = I + 1, P 


IX = 1X1 + J 
X(IX) = 0 


V(2*I) 4 V(2*J-1) 


X(IX) = 1 


GEN 


GRAND=2 


FOR I = 1, P 


RANDU (RAND) 


PCI) = [L3*RAND+1] 


FOR I = 1 to P 


FOR RO = 1 to RI 


FOR DO = 1 to D(I) 


RANDU (RAND) 


R(DO+K(I),R0) = fN3*RND+l] 


CPL 


FOR 

13 = 1 to P 



V(13) = D( I3) . J 

FOR 

I = 1 to P - 1 



1X1 = (I-l)*P-I/2) - I 



FOR J = I + 1 to P 




IX = 1X1 + J 



K X(IX) = 0 




V(J) = MAX(V(J), V(I) 

B38RHI 

KPATH = V(P) 
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1 . 


2 . 


3. 


(INI) Read option 
variables P7,SCH. Set the 
matrix L = matrix A. Set 
vector X to zero. 

(PRNO) Compute the priority 
number by CR1. Store it in 
vector E. 

(COMPET) Get a set of 
simultaneous activities. 

Save the candidates in 
vector S. 

X = 0 if the activity has 
not been scheduled 
X = T1 the time the activity 
is scheduled 
LL = 0 No event between 

activities PO and J 
LL = 1 Otherwise 



M = M + 1 
S(M) => PO 
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4. (ORDER) Sort the set S 
according to their priority 
number, 

5. (SCH) The beginning of 
scheduling the set S.W- will 
keep track of those activities 
that are scheduled or cannot 
be scheduled because of 
precedence constraints. 

T1 is the time counter. 

6. (PRECON) To check if there 

is any precedence constraints. 
IF there is none THEN check 
for resource constraints. 

ELSE go to CHEK. 

7. (ADQ) To check if the 
resources requested by 
activity S(W) is available. 

IF yes THEN schedule that 
activity, and update the 
resource available. ELSE 
go to CHEK. 


ORDER _ 

FOR I = 1 to M - 1 _ 

FOR J = 1 to M ■ I 
K E(SCD) > E(S(I+J)) 
N = SCI) 

S(I) = S(I+J) 


---.——-1 

SCH 


w = 

1 

y\ ti = 

i A 

^ ADQ 

PRECON 


BL = 0 , B2 = 

0 , PO = S(W) 


FOR J = 1, PO - 1 _ 

PREC (LL) 

Nv_ LL = 0 _ 

_ B1 = B1 + 1 

\ X(J) = 0 

X(J) + D(J) > T1 

_ 1 B2 = B2 + 1 

_ B2 = B1 

CHEK 


FOR DO = 1 to D(S(W)) _ 

FOR RO = 1 to RI _ 

YX R(DO+K(SW),R0 ) <_L(T1+D0-1,R0) /N 

_ 1 -*• CHEK _ 

X(SW) = T1 , V(J2) = S(W) , J2 = J2 * 1 

FOR DO = 1, D(S(W)) _ 

FOR RO = 1, RI _ 

L(Tl+DO-l,R0) = L(T1+D0-1,R0)- 
_ _ R(DCH-K,RO) 
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8. (CHEK) To check if all 
the candidate activities 
are checked. IF_ they are, 
THEN check if any activity 
has not been scheduled. 

IF there are any, THEN 
order the already scheduled 
activities and prepare to go 
to GOBACK. ELSE go to 
OUTPUT. ELSE go to PRECON. 

9. (GOBACK) To prepare the 
variables for going back 
either to schedule whatever 
activities are left in S or 
to find a new competing set 
of activities. Update the 
time clock. 


CHEK 


5 




W = M 

/n 


W = W + 1 


+ PRECON 



FOR PO = 2, P 




Y\ 


X 

*0 

O 

V_ S 

II 

o 

_ A 




FLAG = 1 




FLAG = 1 

_ A 


OUTPUT 

K 




J2 = 0 

A 



FOR I = 1 to J2 - 1 





FOR J = I + 1 to J2 - 1 






vm > v(j) 

_ A 





N = V(I) 

v(i) = v(j) 

V(J) =■ N 


GOBACK 


CO = 1 


FOR I = 1 to M 


Y'V X(S(I))> o A 


Cl = Cl + 1 

\ Cl 4 m A 


T1 = T1 + 1 


SCH 

K3 = 0 


FOR I = 1 to M 


y\ x(S(i)) = o A 


K3 = K3 + 1 

Y\ K3 4 1 A 

T1 = MAX(T1,X(S(I)) T1 - T1 + 

♦ D(S(D) D(S(D) 

COMPET 











(OUTPUT) To print the 
calculated schedule. 

PREC - A block used to 
calculate the vector image 
of the adjacency matrix. 
Then pass LL. 

MAIN - The main program 
that calls the subroutines 
SETUP and HSP. 


OUTPUT 



1X1 = (J-l)*(P-J/2) - J * PO 

_ LL = XB(IXI) _ 

RETURN 


SETUP 
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